The N = 16 spherical shell closure in 24 
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The unbound excited states of the neutron drip-line isotope 24 have been investigated via the 
24 0(p,p') 23 0+n reaction in inverse kinematics at a beam energy of 62 MeV/nucleon. The decay 
energy spectrum of 24 0* was reconstructed from the momenta of 23 and the neutron. The spin- 
parity of the first excited state, observed at _E X = 4.65 ± 0.14 MeV, was determined to be J w = 
2 + from the angular distribution of the cross section. Higher-lying states were also observed. The 
quadrupole transition parameter P2 of the 2^ state was deduced, for the first time, to be 0.15 ± 
0.04. The relatively high excitation energy and small P2 value are indicative of the N = 16 shell 
closure in 24 O. 



PACS numbers: 21.10.Re, 25.40.Ep, 27.30.+t 

Magic numbers are a unique feature of finite Fermionic 
quantum systems. In the case of atomic nuclei, exper- 
iments over the last decade or so have shown that the 
well-known magic numbers (2, 8, 20, 28, . . . ) seen to 
occur in stable nuclei often disappear and are replaced 
by new ones as the neutron or proton drip lines are ap- 
proached [1]. For instance, the low excitation energies 
(E x ) and large quadrupole transition parameters (fa) 
of the first 2+ states in ^ 2 Be 8 [2] and ?2 M §20 [3] point 
to the disappearance of the magic numbers N = 8 and 
20. Higher in mass, the recent observation of a low-lying 
2^ state in 42 Si28 provided evidence for the "melting" of 
the N = 28 shell closure [4]. 

In the oxygen isotopes, new shell closures at N = 14 [5- 
7] and 16 [8-16] have been proposed. A high excitation 
energy and low B(E2) value of the 2+ state of 22 O was 
determined via Coulomb excitation [5] , which is sensitive 
to the charge distribution. Proton inelastic scattering on 
22 O, which is sensitive to both the proton and neutron 
distributions, was studied by Becheva et al. [7] and a 
small quadrupole transition parameter /? 2 (0+ s — >2+) = 
0.26 ± 0.04 was reported. Taken together these results 
indicate the near spherical character of 22 O and a sizable 
gap at N = 14. Recently, the low-lying level structure of 
24 was studied via proton knock-out from 26 F by Hoff- 
man et al. [14]. They reported an even higher excitation 
energy (E x = 4.72 ± 0.11 MeV) for the first excited state 
of 24 than that of 22 (3.20 ± 0.01 MeV) [6] suggesting, 



as proposed by Brown and Richter [16], a shell closure at 
N = 16. It should be noted, however, that the spin-parity 
assignments for the 24 excited states were only based 
on a comparison with the predicted shell-model energies. 
Otsuka et al. have investigated theoretically the struc- 
tural evolution of the oxygen isotopes with increasing 
neutron number (iV) and attributed the development of 
the shell closure at N = 16 to the strong neutron-proton 
tensor interaction [10, 11]. 

In this Letter we report on the first spectroscopic study 
of 24 O by proton inelastic scattering. In addition to the 
excitation energies of the states populated, the well- 
known character of proton inelastic scattering also per- 
mits the spins-parities, as well as the quadrupole transi- 
tion parameter (fa) of the first 2+ state to be deduced. As 
described below, we have been able to provide a firm 2 + 
assignment for the state at E x — 4.65 ± 0.14 MeV and 
determine the fa > the small value of which is indicative of 
the spherical closed-shell character of 24 O. A comparison 
of the E x (2j) and fa for the chain of oxygen isotopes 
shows strong evidence for a large shell gap at N — 16. 

The experiment was performed at the RIPS facil- 
ity [17] at RIKEN. A schematic view of the down- 
stream section of RIPS and the experimental setup is 
shown in Fig. 1. The 2 O secondary beam was produced 
using a 1.5 mm-thick Be production target and a 95 
MeV/nucleon 40 Ar primary beam with a typical inten- 
sity of 40 pnA. The intensity of the 24 O secondary beam 
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was ^4 ions/sec with a momentum spread of Ap/p = 
± 3%. The liquid-hydrogen (LH 2 ) target [18] was in- 
stalled at the achromatic focus F3 of RIPS. The effective 
target thickness and the mid-target energy were 159 ± 
3 mg/cm 2 and 62 MeV/nucleon, respectively. The sec- 
ondary beam was tracked particle-by-particle on to the 
target using two drift chambers (NDCs) located just up- 
stream of the target. The Bp-TOF-AE method was em- 
ployed to identify the charged fragments following reac- 
tions of the 24 beam with the LH 2 target. The magnetic 
rigidity (Bp) was determined from the position and an- 
gle measurements derived from the MDC and FDC drift 
chambers (Fig. 1). The TOF (time-of-flight) of the frag- 
ments was measured between the target and the plastic 
scintillator charged particle hodoscope (CC). The energy 
loss (Ai?) was measured using the CC. 

The neutrons were detected using the plastic scintilla- 
tor neutron counter (NC) array placed some 4.7 m down- 
stream of the target together with the charged particle 
veto counter (Fig. 1). The NC array consisted of four 
layers with a total thickness of 24.4 cm. The momentum 
vectors of the neutrons were determined from the TOF 
between the target and the NC along with the hit posi- 
tions at the target and the neutron counter. A neutron 
detection efficiency of 25.0 ± 0.8% at 64 MeV was mea- 
sured for a 2 McVec threshold in a separate 7 Li(p, n) run. 

The decay energy spectrum of 24 0* was reconstructed 
from the measured four momenta of 23 O and the emitted 
neutron. The decay energy, -Edccay, is expressed as: 

^dccay = \J {Ef + E n ) 2 - \pj + p n \ 2 — (M f + M n ) , (1) 

where Ef(E n ) and Pf(p n ) are the total energy and the 
momentum of 23 O (neutron) and Mf and M n are the 
masses of 23 and the neutron, respectively. 

Figure 2(a) shows the coincidence yield for 23 O and a 
neutron in decay energy. Two peaks are clearly visible at 
-^dccay ~ 0.7 and ~3.2 MeV. Figure 2(b) shows the decay 
energy spectrum in terms of cross section (dcr/d-Edecay) 
after correction for the detection efficiencies and accep- 
tances. The error bars are statistical only. The geometri- 
cal acceptance was estimated using a Monte Carlo simu- 
lation taking into account the beam profile, the geometry 
of the setup, the experimental resolutions, and the mul- 
tiple scattering of the charged particles. The acceptance 
of the neutron detector drops rapidly as the decay en- 
ergy increases, resulting in a significant suppression of 
the higher-lying strength. 

The decay energy spectrum was fitted using two res- 
onance distributions for the first peak at -Ed 0C ay ~ 0-7 
MeV, a Gaussian distribution for the broad feature at 
-^decay ~ 3.2 MeV, and a Maxwell distribution for the 
non-resonant continuum background [19]. The first peak 
around ~0.7 MeV was interpreted as two closely spaced 
resonances because the width is greater than the es- 
timated experimental energy resolution by a factor of 




FIG. 1. A schematic view of the experimental setup. 
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FIG. 2. (Color online) (a) The coincidence yield of 23 
and a neutron in decay energy and (b) the cross section, 
dcr/diSdccay. The error bars are statistical only. The red-solid 
and blue-dashed histograms are the results of fits for reso- 
nances at -Bdecay = 0.56 MeV and 1.06 MeV, respectively. The 
dotted Gaussian function at E^ CCBjy ~ 3.2 MeV represents 
higher-lying state(s). The dot-dashed line corresponds to the 
non-resonant continuum [19]. 

~2. The experimental energy resolution is energy de- 
pendent and estimated to be AEdccay ~ 0.5 a/ '-Edccay 
(MeV) in FWHM. The single-particle width for a d- 
wave neutron resonance is predicted to be less than ~0.1 
MeV [20]. Therefore, the experimental resolutions were 
adopted for the widths. The solid histogram of Fig. 2(b) 
represents the best fit result of the overall distribution of 
the total cross section, dcr/dE'decay, which consists of two 
resonances of the first peak, a Gaussian distribution for 
the higher lying peak at ~3.7 MeV and the non-resonant 
continuum background. The red-solid and blue-dashed 
histograms are those for the first two resonance states 
at £ d ccay = 0.56 ± 0.05 MeV and 1.06 ± 0.10 MeV, 
which correspond to E x = 4.65 ± 0.14 MeV and 5.15 ± 
0.16 MeV, respectively, adopting the separation energy 
S„ = 4.09 ± 0.13 MeV [21]. The dot dashed line is the 
distribution for the non-resonant continuum background 
obtained from the same fit. The excitation energy of the 
first resonance is consistent with that of the recent study 
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FIG. 3. (Color online) Angular distributions for the reso- 
nances at -Bdccay = 0.56 (filled circles) and 1.06 MeV (open 
circles in the inset). The error bars are statistical only. The 
curves represent the results of microscopic (black curves) and 
phenomenological (red curves) DWBA calculations (see text). 
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FIG. 4. (Color online) The neutron number dependences of 
(a) E*(2f) and (b) fa (2f)- The present results are shown 
by the filled circles. The -E x (2+) for 7V= 10-14 were taken 
from Refs. [6, 34]. The fa for N= 10-14 were taken from 
Refs. [7, 35]. The dotted lines represent the USDB shell model 
predictions [24]. 



of Hoffman et al. (4.72 ± 0.11 MeV) [14] and in accord 
with the vlsx/2 - v0d 3/2 shell gap (4.86 ± 0.13 MeV) 
derived from the location of the 25 O g . s . resonance [13]. 

The dotted line in Fig. 2(b) shows the fit to the high- 
lying peak with -Edecay ~ 3.2 MeV corresponding to an 
excitation energy of 7.3 MeV. Recently a high-lying state 
has been reported at E x ss 7.5 MeV, which was produced 
via proton knock-out from 26 F and de-excited to 22 O g . s . 
via two neutron cascade, 22 O g . s .+2n [22]. This is a differ- 
ent decay channel from the present one, 23 O g . s .+n. Since 
the decay channel is different, the strength observed here 
is additional to that previously reported, thereby imply- 
ing considerable strength in the excitation energy region 
around 7.3 MeV. The nature of this high-lying strength 
will be discussed later. 

The total inelastic cross sections to the low-lying reso- 
nances at .Edocay = 0.56 and 1.06 MeV were determined 
to be a = 2.6 ± 1.1 mb and 2.2 ± 1.2 mb, respec- 
tively, after subtracting the non-resonant continuum. The 
quoted errors mainly come from the uncertainties in the 



fitting (—39% and —47%) and the choice of the functional 
form describing the non-resonant continuum (—19% and 
-27%). 

The single-particle configurations for the first and sec- 
ond excited states of 24 are expected to be mainly 
v{ls 1/2 )~ 1 v{M ?i /2) 1 J* = 1+ or 2+ in the shell model 
picture. Calculations using the universal sd-shell model 
interactions, USD [23] and USDA/B [24], as well as that 
including the continuum states [25] all predict that the 2 + 
state is lower lying than the 1 + state by approximately 
0.5-1.0 MeV. 

The angular distribution, dcr/df2, for the first reso- 
nance state of Fig. 3 was obtained by fitting the decay 
energy spectrum of each angular bin (of 20°) in the simi- 
lar way to the fitting method explained for Fig. 2(b) but 
adopting the uniform angular distribution for the non- 
resonant continuum background. With this the decay en- 
ergy spectrum at each angle could be fitted with only 
the resonance distributions after the subtraction of the 
non-resonant background whose total cross section was 
fixed in the explained procedure for Fig. 2(b). The error 
bars are statistical only. The experimentally determined 
angular distributions have been compared with micro- 
scopic DWBA calculations performed using the dw81 
code [26]. The global optical potential KD02 [27] and 
the M3Y [28] were employed for the distorted wave func- 
tion and for the effective nucleon-nucleon interaction, re- 
spectively. The one-body transition densities were cal- 
culated using the shell- model code nushell [29]. The 
calculations used the USDB interaction for the 2^ and 
the 1^ states, and the WBT [30] interaction for the 3± 
state. The size parameter b — 2.03 ± 0.08 fm was chosen 
to reproduce the rms radius of 24 [31] within the har- 
monic oscillator potential. The results of the calculation 
for the transitions g s — >2+ (black-solid line), 0+ s — >lf 
(black-dotted line), and + s — (black-dashed line) are 
shown in Fig. 3. It may be noted that there are no ad- 
justable normalization parameters in the calculation. The 
result for the + s — ^2^ {v{lsi/ 2 )~ 1 i'(0d 3 /2) 1 } transition 
reproduces the angular distribution of the resonance at 
-^decay = 0.56 MeV very well, strongly supporting a spin- 
parity assignment = 2 + . 

The quadrupole transition parameter /3% is a measure 
of the deformation [7]. We have derived the f3 2 value of 
the 2^ state of 24 O by normalizing the phenomenolog- 
ical collective DWBA calculation to the measured total 
cross section. The red-dashed and red-dot-dashed lines 
in Fig. 3 represent the normalized angular distributions 
of the collective DWBA calculations using the ECIS97 
code [32] for the two different optical potentials KD02 
and CH89 [33], respectively, from which we deduce f} 2 
= 0.15 ± 0.04. The error reflects the uncertainty in the 
total cross section and the choice of optical parameters. 

Figure 4 shows the neutron number dependences of 
the measured and calculated (dotted lines) E x (2f) and 
/?2(2^ ) for the even-even oxygen isotopes. The E x (2f) 
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values increase considerably in moving from N = 12 to 
16. This trend and the high excitation energy at N = 
16 are well reproduced by the USDB shell model calcula- 
tions. In the case of 24 (N = 16), the 2^ state lies ~1.5 
MeV higher than in 22 0, reflecting the large gap between 
the ulsi/2 and vOd 3 / 2 orbitals, which is a maximum at 
N = 16. 

The dotted line in Fig. 4(b) represents the USDB shell 
model predictions using effective charges of e p = 1.36 
and e„ = 0.45 [36] and following Bernstein's prescrip- 
tion [37] . The calculations reproduce reasonably well the 
experimental j3 2 values, in particular those of 22 O and 
24 0. Based on the small f3 2 value (0.26 ± 0.04) de- 
rived from the proton inelastic scattering of 22 O, Becheva 
et al. concluded that the ^0g? 5 / 2 sub-shell is closed at N 
= 14. For 24 0, the f3 2 value (0.15 ± 0.04) is even smaller 
than that of 22 O, thereby suggesting that it is the least 
deformed of the oxygen isotopes. From the E x (2~£) and 
(3 2 values of 22 O and 24 O, it is clear that the gap at N 
= 16 must be considerably larger than at TV = 14. This 
is in accord with the predictions of Otsuka et al. [10-12] 
whereby the vls^n sub-shell is closed in 24 0. 

While the cross section for the first excited state is well 
reproduced by the microscopic DWBA calculations for 
the Q+ s - >2j {v(lsi/ 2 )~ 1 v(0d 3 / 2 ) 1 } transition, that of 
the second resonance, 2.2 ± 1.2 mb, is much larger than 
that calculated for the 0+ s — sT^ {v(ls 1 / 2 )~ 1 u(0d3/2) 1 } 
transition (Fig. 3). While the uncertainty in the measured 
cross section (56%) is large, one conjecture is that some 
of the missing strength might be attributed to negative- 
parity states v{lsi/2)~ 1 v{fj>) 1 that come down in energy 
to lie near the 2+ state owing to the quenching of the gap 
between the v0d 3 / 2 orbital and the vfp shell. Indeed, low- 
lying intruder 3/2]~ and 7/2]~states in 27 Ne [38] indicate 
a narrowing of the gap between the vQd 3 / 2 orbital and 
the vfp shell. 

As mentioned earlier, the high-lying strength (E x ~ 
7.3 MeV) observed here represents an additional contri- 
bution to that observed previously in the two-neutron 
emission channel [22] . High- lying states produced via the 
promotion of a 0d 5 / 2 neutron would decay most likely 
to the 5/2 + first excited state of 23 0— a v{0d 5 / 2 )^ 1 hole 
state— which is unbound [39, 40] and, in turn, decays to 
22 0g.s. via neutron emission, a process not observable 
in the present experiment. This implies that the states 
observed here around ^7.3 MeV would have mostly 
v{^si/ 2 )^ 1 v{fp) 1 negative-parity configurations. 

In summary, we have investigated the unbound excited 
states of neutron-rich 24 via proton inelastic scatter- 
ing in inverse kinematics. The excitation energy of the 
first excited state was determined to be 4.65 ± 0.14 MeV 
and the spin-parity was assigned J x = 2+. In addition, 
a relatively small j3 2 parameter was determined (0.15 ± 
0.04), indicative of the spherical character of 24 and 
the large shell gap at N — 16, confirming theoretical 



predictions [10-12, 16]. Finally, the strong higher-lying 
strength, identified as negative-parity v{^si/ 2 )^ l v{fp) 1 
excitations, is suggestive of the quenching of the gap 
between the neutron sd and fp shells, as observed in 
27 Ne [38]. 
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